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P^ I Abstract. In the presence of large-scale vortical motions and/or magnetic-field strains, the 

^0 ' turbulent cross helicity (velocity-magnetic-field correlation in fluctuations) may contribute to 

(-H ' the turbulent electromotive force and the Reynolds stress. These effects of cross helicity are 

O „ considered to balance the primary effects of turbulence such as the turbulent magnetic diffusivity 

I ■ in magnetic-field evolution and the eddy viscosity in the momentum transport. The cross- 

O ' helicity effects may suppress the enhanced transports due to turbulence. Physical interpretation 

of the effects is presented with special emphasis on the difference between the cross-helicity effect 
and the usual a or helicity effect in the dynamo action. The relative importance of the cross- 
helicity effect in dynamo action is validated with the aid of a direct numerical simulation (DNS) 
of the Kolmogorov flow with an imposed magnetic field. Several mechanisms that provide 
turbulence with the cross helicity are also discussed. 



1. Introduction 



Cross helicity (velocity-magnetic-field correlation) is a quantity of fundamental importance at 
high magnetic Reynolds number magnetohydrodynamic (MHD) flow. The total cross helicity, 
as well as the MHD energy, is an inviscid invariant of the MHD system of equations, but a 
pseudoscalar which represents breakage of mirror-symmetry. Cros s helicity is expected to play 
K> I a certain role in the turbulent dynamo (see lYoshizawa et all 12004 . and works cited therein) . 



j^ ' In the presence of a large-scale vortical motion, the turbulent cross helicity {u' ■ h') contributes 

to the electromotive force aligned with the large-scale vorticity fl{= V x U) {u': velocity 
fluctuation, b': magnetic-field fluctuation, U: mean velocity, (•••): ensemble average). A 
fluid element in the large-scale vorticity is subject to the Coriolis-like force u' x f2. Provided 
that there is a positive (or negative) cross correlation between the velocity and magnetic-field 
fluctuations, a contribution to the electromotive force parallel (or antiparallel) to the mean 
vorticity, {5u' x b') ex {u' ■ b')fl, arises (figured]). This is in marked contrast with a positive (or 
negative) turbulent helicity {u' ■ u'), correlation between the velocity and vorticity fiuctuations 
may contribute to the electromotive force antiparallel (or parallel) to the mean magnetic field 
[a;'(= V X u'): vorticity fluctuation]. 
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If we take the cross- helicity effects into account, the turbulent electromotive force -Em is 
expressed as ( Yoshizawal . ll99Cl ) 



Em = iu' X b') = aB - /3V X B + 7/2, 



(1) 



where a, /3, and 7 are the transport coefficients expressed in terms of the residual helicity 
{-u' -u' + b' ■ j'){= H), the turbulent MHD energy (w'^ + b'^)/2{= K), and the cross helicity 

{u' ■ b'){= W), respectively. 
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Figure 1. Physical interpretation of the cross-helicity effect. Redrawn from lYokoil ( 19991 ). 



If we substitute ([T]) into the mean magnetic induction equation, we obtain 

-_ = Vx(C/xS) + Vx (aB + 7/?) - V X [(?? + /3)V x B] 
ot 



(2) 



(rj: magnetic diffusivity). First, in the presence of turbulence, the effective magnetic diffusivity 
or resistivity is enhanced by turbulence (77 — >■ r; + /3) with spatiotemporal variation of /3. So, the 
coefficient /? is called the turbulent magnetic diffusivity or anomalous resistivity. Secondly, the a- 
and 7-related terms express possible magnetic-field generation mechanisms due to pseudoscalars. 
If turbulence has some asymmetric properties represented by these pseudoscalars, the enhanced 
transport due to /3 can be suppressed by these pseudoscalar effects. This suppression of turbulent 
transport is based on the dynamic balance between two turbulent effects: enhancement and 
suppression of the transport. This state is entirely different from the one subject only to the 
molecular transport. Turbulent itself is very strong, but due to the additional asymmetry 
represented by pseudoscalars, effective turbulent transport is dynamically suppressed. Then 
important point to see is how and how much such pseudoscalars can be present in turbulence. 
In the long history of the dynamo study, special attention has been paid to the a or helicity 
effect [the first term in ([T|)]. In contrast, the 7 or cross-helicity effect [the third term in ([1]) 



has been almost missing except for some limited number of works (e.g., lYoshizawa et al. 



and works cited therein). One of the reasons of this missing may be attributed to the current 
treatment of turbulence; in most works, the large-scale velocity is neglected in the consideration 
of turbulence because of the Galilean invariance of the momentum equation. However, such 
treatment inevitably leads to the neglect of the large-scale shear effects. On the contrary, a large- 
scale rotation is ubiquitous and sometimes essential ingredient in astro/geophysical phenomena. 
In this sense, we should take effects of large-scale inhomogeneous flows into consideration of the 
realistic turbulent dynamo studies. 

The number of works that c ontain numerical simulations exaniining or validating the cross- 
helicity effe c t is very limited (iHambal. Il992l : iNishino fc Yokoi Il998l : lYokoi fc Hambal . l2007l : 
Yokoi et all 120081 : ISur fc Brandenburd . l2009l ). In the present work, we examine the effect 



of turbulent cross helicity, and see the validity of the notion of cross-helicity dynamo. This 
must be a very interesting contribution to the astro/geophysical dynamo studies from the shear 
turbulence study. 

Organization of this paper is as follows. After an outlined presentation of the turbulent 
cross-helicity generation mechanisms with the aid of the evolution equation of the turbulent 
cross helicity in section [21 a brief description of the MHD Kolmogorov flow is given in section [3j 
Some of the basic numerical results are presented in section HI which include the analysis of 
the turbulent electromotive force (section I4.ip , the comparison of the induced fields with some 
approximate dynamo solutions (section I4.2p , the comparison of the spatial distributions of the 
turbulent cross helicity and its production rate (section I4.3p , and the turbulent cross helicity 
scaled by the turbulent MHD energy (section 14. 4p . Concluding remarks with directions in the 
future work are presented in section [SI 

2. Cross-helicity generation mechanisms 

From equations of the velocity and magnetic-field fluctuations, u' and 6', we obtain the equation 
for the turbulent cross helicity 

W = {u' ■ b') (3) 

as 

DW ( d \ 

= \- + U ■V\w = Pw-ew^Tw. (4) 



W = -7e'^'7r3r-^M-/?, (5) 



Dt \dt 

Here, Pw, ew, and T\y are the production, dissipation, and transport rates of W. They are 
defined by 

, dB ' 

, /du"'db'^\ 
^^ = (" + ^)W5^/' (') 

Tw = B-VK + V-T'w, (7) 

where V • T^ is the transport rate of W that can be written in the divergence form. It comes 
from the higher order terms of u' and b', whose details are suppressed here. In ([5]), 7?. and 
-Em ai'e the Reynolds stress (combination of the Reynolds and turbulent Maxwell stresses in the 
MHD case) and the turbulent electromotive force, respectively. They are defined by 

and ([1]). The Reynolds stress is known to be expressed as 

7^"/5 = ^Kn6''^ - z/k5"^ + mM''^, (9) 

where Kji{= {u'"^ — 6'^)/2) is the turbulent MHD residual energy, and ^'k[= (7/5)/3] and 
um[= (7/5)7] are the turbulent transport coefficients. The strain rates of the mean velocity 
and magnetic field, S and A^, are defined by 

5"/^ = ^ + ^ (10) 

dxP ^ dx» ■ ^ ' 



If we substitute the expressions for TZ and Em, © and ([T]), into the production rate of the 
turbulent cross hehcity, Pw, dS]), we obtain 



1 1 2 

Pw = -j^kS : M. - -vmM. 



aB ■ n + /3J ■ n --ff2^ 



(12) 



[S : A4 = S"'''M"''', A4^ = M'^^M'^'']. This shows that combinations of the mean-field shears 
coupled with the turbulent transport coefficients give the generation mechanisms of the turbulent 
cross helicity. For the detailed meaning of each term in the production, dissipation, and transport 
rates, and for the possible si tuations where th e cross helicity is supplied to turbulence, the reader 
is referred to lYokoil (120111 1: lYokoi k Hoshind (|201lh . 



3. Magnetohydrodynamic Kolmogorov flow 

Kolmogorov flow is a three-dimensional periodic flow with external forcing. This flow is 
homogeneous in the x and z directions, but due to the external forcing expressed as 



/ = (r , /^ n = (/o sm (27ry/L^) , 0, 0) 



(13) 



(L^, L^, L^: box dimension), it is inhomogeneous in the y direction (figure [2]). The Kolmogorov 
flow provides a good test for investigating three-dimensional inhom ogeneous tu r bulent flow 



simultaneously with the mean-velocity shear and the anisotropy (|Sarris et all 120071 . and 
works cited therein, also see Hamba, 1992). In order to examine the basic properties of 
magnetohydrodynamic (MHD) turbulence, in addition to the external forcing, we impose a 
uniform large-scale magnetic field in the inhomogeneous (y) direction as 



B 



(i?^s^s^) = (o,5o,o). 



(14) 




A f(y)=fo sin (2ny/Ly) 
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Figure 2. Kolmogorov flow with a magnetic field imposed. 



In this work, with the aid of a direct numerical simulation of the Kolmogorov flow, we address 
the question of whether the cross-helicity effect is relevant to the turbulent electromotive force 
-Em or not. 

It was reported that the statistics of turbu lence depend on the several computational 
conditions in the Kolmogorov ( Sarris et a/.l . 120071 ) . In this sense, to obtain quantitative results. 



we have to scrutinize these conditions, including the dependences on the Reynolds number, 
aspect ratio of the box, strength of the imposed magnetic field, etc. We also have to take care 
of the characteristic length scales in the three directions (in particular in the direction of the 
imposed magnetic field) in comparison with the domain sizes. However, here our attention 
is focused on the basic behaviours of the cross- helicity effect in Em- So, we only consider a 
situation of a cubic box with the numerical discretization of 256^. In the following, the operator 
(■ ■ ■ ) denotes averaging over the homogeneous (x and z) directions and ensemble averaging over 
70 independent realizations in time. The primed quantities correspond to deviations from the 
averages. Each realization shows large spatial variations. In order to get more smooth statistics, 
we have to take larger statistical ensembles. 

4. Basic results and discussions 

The situation of Kolmogorov flow may be unrealistic in that it is difficult to experimentally 
generate such a periodic forcing (|13p in an unbounded flow. Nevertheless, this flow gives a 
test for investigating three-dimensional turbulence with inhomogeneity, anisotropy, and velocity 
shear. So, with an imposed magnetic field, it gives a suitable situation to validate the cross- 
helicity effect in a magnetohydrodynamic turbulent flow. In this work, we examine mainly two 
points: (i) how the cross-helicity-related term contribute to the turbulent electromotive force; 
and (ii) how the cross helicity is supplied to turbulence. 

4.I. Turbulent electromotive force 

First we consider the turbulent electromotive force Em_{= {u' x b')). The spatial distribution of 
the z component of Em is shown in figure [3l The profile of Em is sinusoidal as expected from 
the external forcing. 
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Figure 3. Spatial distribution of the turbulent electromotive force and each dynamo term in 
the Kolmogorov flow versus y: {u' x b')^; aB^; l3J^', - ■ - 7^7^. 



In figure [3l we also showed the spatial distribution of each term of the right-hand side in ([T|) . 
We see from this figure that the a- or helicity-related term is much less than the /3- or turbulent 
magnetic diffusivity-related term and 7- or cross-helicity-related terms. The main contribution 



to the turbulent electromotive force in this case is attained by the balance between the /3- and 
7-related terms: 

Em^-^J + I^. (15) 

Namely, the turbulent magnetic-diffusivity effect is balanced by the turbulent cross-helicity effect 
rather than by the turbulent helicity effect. 

This result certainly shows the importance of the cross-helicity effect as compared to the a 
or helicity effect in the dynamo action in a shear flow situation. 

4-2. Dynamo solutions 

If the a or helicity effect is negligibly small in the turbulent electromotive force, the mean 

induction equation ([2]) is written as 

= Vx [C/xS-(/3 + r/)J + 7/2]. (16) 

ot 

In the present geometry of the Kolmogorov flow, this reads 

QQX QJJX Q Q 

The first term in ()17p arises from the velocity inhomogeneity in the direction of the mean 
magnetic field, which corresponds to the so-called effect. 

Fo r a stationary state , we may consider an approxi mate dynamo solution for the induced 
field (JYoshizawa fc Yokoil . Il99.i lYokoi fc Hoshinol . 120111 ) as 



7 W 

SB = jU = C-U, (18) 

where C is the model constant of 0(10^^). Here we dropped the magnetic diffusivity rj as 
compared with the turbulent counterpart /3. However, in case 7//? or W/K spatially varies, the 
approximate dynamo solution (fT8]) is not the case. In figure [D^a) we plot the induced magnetic 
field B^ in comparison with the approximate dynamo solution {W / K)U^ . 

If the turbulent magnetic diffusivity is fully balanced by the cross-helicity effect in the 
turbulent electromotive force, the mean electric-current density is aligned with the mean vorticity 
as 

§j = l{2 = C—n. (19) 

p K 

We plot the induced mean electric-current density J^(= —dB^/dy) in comparison with the 
approximate dynamo solution {W / K)Q.^[= —{W/K){dU^/dy)\ in figure [D^b). Unlike (fT8|) . the 
general tendency of J^ is in agreement to some extent with the counterpart of {W/K)Vt'^, 
although a further analysis of each term of (|17p is needed. 

4-3. Cross helicity and its production rate 

Next, we consider the cross-helicity generation mechanisms in the Kolmogorov flow. For this 
purpose, we show in figure [5] that the spatial distribution of the turbulent cross helicity scaled by 
the characteristic time of turbulence, W/t. As for the time scale, we adopt the energy cascade 
time defined by the turbulent MHD energy K divided by its dissipation rate e as 

r = K/e (20) 

with 

K = (-u'2 + 5'2^/2, (21) 
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Figure 4. Spatial distribution of the induced field and approximate dynamo solutions versus 

y. (a) Magnetic field S^ and velocity {W/K)U'': (black) fi^; (gray) {W/K)U''. 

(b) Electric-current density J^ and vorticity {W/K)VL'^: (black) J^; (gray) (W/K)VL'^ 

versus y. Constant factor (model constant) C is not included in the dynamo-solution expressions. 
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We also plot each term in the production rate of the turbulent cross helicity, 
figure O 



(22) 
and d?]), in 
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Figure 5. Spatial distribution of the turbulent cross helicity and its production-rate terms in 
the Kolmogorov flow versus y: W/t; -W^idB" /dx^); ^m • ^; B VK. 



This figure shows that the spatial distribution of the turbulent cross helicity is in good 
agreement with the counterparts of the cross-helicity generation mechanisms. Of course we 



need to delve further into the detailed budget of the cross-helicity evolution dH) including the 
dissipation and transport rates. With this reservation in mind, we see from this result that the 
production rate provides a good measure for the turbulent cross helicity itself. 

Estimate of the turbulent cross helicity with its spatial distribution from observation is a 
challenging problem. All the three components of the velocity and magnetic-field fluctuations 
should be measured simultaneously. In actual observations, we often encounter a situation where 
the direct measurement of the cross helicity is impossible or hard to achieve. Even in such a 
case, if we can measure the production rate of cross helicity, we may estimate the turbulent cross 
helicity. In the detailed expression (|12p for the production rate of W, the second and fifth terms, 
—umM. /2 and —-yfi , always work for the reduction of W. If we concentrate our attention on 
the genuine generation mechanisms of W, which come from the primary effects of turbulence 
such as the turbulent-energy- or /3- and i^K-related terms, then we have 



Pw^^i^kS -.M + ^J ■ O. 



(23) 



In this case, without knowing any information on the turbulent cross helicity itself, we may be 
able to estimate W by measuring the large-scale field configuration represented by the mean 
velocity and magnetic-field strains, S and A^, the mean electric current density J, and the 
mean vorticy f2, etc. 

4-4- Scaled cross helicity 

As we see in (|19p . the proportional coefficient j/f3 or W/K (the turbulent cross helicity {u' ■ b') 
scaled by the turbulent MHD energy {u''^ + &'^)/2) combines the large-scale vorticity /2 with 
the large-scale electric-current density J induced by the dynamo action. If the turbulent cross 
helicity is negligibly small compared to the turbulent MHD energy, we have substantially no 
dynamo effects due to the cross helicity. In this sense, W/K is one of the most important 
quantities in the arguments of the cross-helicity effects. 

The spatial distribution of the scaled cross helicity W/K versus y is plotted in figure El 
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Figure 6. Spatial distribution of the turbulent cross helicity W{= {u' ■ b')) scaled by the 
turbulent MHD energy K{= (ti'^ + b'^)/2) in the Kolmogorov flow versus y. 



We see from this figure that the magnitude of the scaled cross helicity is 0(10 ^) — 0(10 ^) 
in this case. It is worth noting here that the absolute value of W/K is mathematically bounded 



as 

K ~ (w'2 + b'2)/2 ^ '■ ^^^^ 

One self-evident but important consequence of this boundedness is that the turbulent cross 
helicity W can not be sustained in the absence of the turbulent MHD energy {K = 0) . 

As the W evolution equation (j3|) shows, the value of W/K must depend on the magnitude 
of the velocity shear and especially on the magnitude of the imposed magnetic field, Bq. To 
obtain more quantitative understanding, we have to delve into the detailed budgets of both the 
turbulent cross-helicity and MHD-energy evolutions. The highest value of \W\/K is realized 
only when the velocity fluctuation u' is totally aligned with the magnetic-field fluctuation b' 
with equipartition of the turbulent kinetic and magnetic energies ( Alfvenic state) : 

u' = ±b' (25) 

This state is often observed in the solar-wind turbulence near the Sun because of the dominant 
effect of the Alfven waves. On the other hand, in the analysis of the galactic dynamos using the 
cross-helicity effect, it was inferred that \W\/K of O( 10~^) — 0(W ~^) is large enough to explain 
the magnetic- field strength of several spiral galaxies ( Yokoi 119961 ). 



5. Concluding remarks 

On the contrary to the current understanding on the turbulent dynamo, where the a effect 
combined with the differential-rotation effect (called the fi effect) plays the main role in the 
generation and evolution of the dynamo magnetic field, the relative importance of the cross- 
helicity effect to the a effect is suggested in the present work. This may pave a new way for the 
turbulent dynamo model for the astro/geophysical context. At least, the cross-helicity effect, in 
other words, the effect of the large-scale shear flow, should be taken into account in the dynamo 
study as a complementary effect to the well-known a or helicity effect. 

A large-scale magnetic field alters statistical properties of turbulence compared to the 
hydrodynamic turbulence. For the sake of further quantitative analysis, the dependences of 
the simulation results on the magnitudes of the imposed magnetic field, velocity shear, and on 
the grid size may be important. These points should be further included in the forthcoming 
paper. 

As a final remark, we refer to the cross-helicity effect in the momentum equation. In the 
momentum equation, the turbulent transport is expressed by the Reynolds stress ([8]). Since 
the cross-helicity or z/M-related term appears in the Reynolds-stress expression ([9]) as a balancer 
with the the eddy viscosity i/k, the cross helicity is expected to play an important role also 
in the suppression of turbulent linear- and angular-momentum transports. In addition, the 
cross-helicity effect may play important role in the momentum transport through the mean 
Lorentz force. The mean electric current J configuration induced by the cross-helicity effect is 
in general not aligned with the large-scale magnetic field (J jf B). This is in marked contrast 
with the counterpart by the a or helicity effect, where J is aligned with B {J x B = 0, 
force-free field). This feature o f the cross - helici ty effect was fully utilized in the studies of 
the inner osci llations in the Sun ( Itoh et a/.l . l2005l ) and flow-turbulence interaction in magnetic 



reconnection ( Yokoi Sz Hoshind . l201ll ). 



In the context of solar physics, the differential rotation of the fluid inside the Sun has 
been left as an unsolved big problem. Recent elaborated numerical simulations showed that 
the Reynolds stress plays a k e y role in the angular-m omentum transport inside the Sun 



tne ileynoias stress plays a K e y roie m tne angular-m omentum transport msiae tne bur 
( Brun. Miesch &: Toomrd . |200J; iMiesch Sz Toomrd . l2009l ). In order to address this problem 



we should examine the Reynolds stress and its model expression ([9]). In addition to the analysis 



of the turbulent electromotive force mentioned in section U we will discuss the Reynolds stress 
expressions in the future work. 
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